use "/mi_data.dta" ///dataset generated with "multiple_imputation.R" 
mi import flong, m(imp) id(country_id year id) imp(e_peinfmor_ln-e_peinfmor_ln_rt)
tabulate year, generate(yr)
drop yr1 
mi tsset country_id year
mi xeq: sort country_id year; gen f_e_peinfmor_ln=F.e_peinfmor_ln; gen l_e_peinfmor_ln=L.e_peinfmor_ln

mi estimate , cmdok: xtreg F.e_peinfmor_ln v2x_EDcomp_thick nocorrupt e_migdppc_ln e_miurbani e_migdpgro e_peinfmor_ln_rt e_peinfmor_ln i.year, fe

mi estimate , cmdok: xtscc f_e_peinfmor_ln v2x_EDcomp_thick nocorrupt e_migdppc_ln e_miurbani e_migdpgro e_peinfmor_ln_rt e_peinfmor_ln yr*, fe
